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A nonlinear kinetic exclusion model is used to study osmosis and pressure driven flows through 
nearly single file pores such as antibiotic channels, aquaporins, zeolites and nanotubules. Two 
possible maxima in the steady state flux as a function of pore-solvent affinity are found. For small 
driving forces, the linear macroscopic osmotic and hydraulic permeabilities Poa and Lp, are defined 
in terms of microscopic kinetic parameters. The dependences of the flux on activation energies, pore 
length and radius, and driving forces are explored and Arrhenius temperature dependences derived. 
Reasonable values for the physical parameters used in the analyses yield transport rates consistent 
with experimental measurements. Experimental consequences and interpretations are examined, 
and a straightforward extension to osmosis through disordered pores is given. 
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time average flux [l/t) 




molecular transit times (S,^^) 
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pore length 
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solvent chemical potential difference 
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hydraulic permeability 
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hydrostatic pressure difference 
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number of pore sections 
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osmotic pressure difference 
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osmotic permeability 
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I. INTRODUCTION 

Transport through microscopic pores is vital to many 
biological and industrial processes |-|]. One phe- 
nomenon involving flow through microscopic channels is 
osmosis. Ubiquitous in biological functions such as cellu- 
lar volume control, osmosis and reverse osmosis are also 
exploited in industrial processes such as solution concen- 
tration, flltration, and catalysis. 

Living cells must transport nutrients, ions and water 
across their lipid membranes. Various biological struc- 
tures and machinery, such as active transporters, cotrans- 
porters, antiporters, and simple channels or pores have 
evolved to perform these tasks Of central impor- 

tance in living organisms is the flow of water. Certain 
cells, such as nephrons, have a permeability to osmot- 
ically driven water transport too high to be explained 
by simple permeation through an oily lipid membrane 
bilayer. In addition to generic antibiotic pores such as 
gramicidin, pores which are apparently water specific 
(aquaporin-CHIP) have been identified By inserting 
aquaporins into their membranes, cells can in principle 



control overall osmotic permeabilities. 

Mechanisms of fluid transport in confined geometries 
are also relevant in industrial applications, especially in 
separation processes. Both natural and manmade mate- 
rials such as zeolites contain many molecular-sized, es- 
sentially single-file channels which can selectively absorb 
fluids. This size specificity can be exploited in separa- 
tion of a mixture of linear and branched chain alkanes, 
where the zeolite acts like a sponge absorbing only the 
desired specie(s) Confining particles in zeolite pores 
can also serve to catalyze reactions . Therefore, there 
has been much research into fluid structure and transport 
within confined pores 0-|l^. 

Although the equilibrium properties of osmotic "pres- 
sure" are understood in terms of macroscopic forces p^ ] 
and statistical mechanics p^ , nonequilibrium osmotic 
transport has been less well categorized. One approach 
uses macroscopic flow equations (such as the parabolic 
Poiseuillean fluid velocity proflle through a pipe p9|); 
these, applied to structures of molecular dimensions lead 
to inaccurate prediction of membrane pore density, or 
pore radius, particularly in biological examples P,pO[. 
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Use of macroscopic parameters such as viscosity, or hy- 
draulic flow, to microscopic transport has led to mod- 
els that do not accurately describe most of the relevant 
microscopic physical processes |2^j2^,^ and many con- 
flicting viewpoints j^]. A recent review of the historical 
controversies and misconceptions is given by Guell and 
Brenner JTzt , who present a clarifying macroscopic de- 
scription of osmosis. Finkelstein also mentions the 
many conceptual disagreements associated with osmosis 
and nonequilibrium thermodynamics, particularly when 
applied to water transport across biomembranes. 

In this paper, we consider molecular ly-sized pore that 
are nearly single-file. Many natural examples of pores 
are of molecular sizes. Zeolites typically have pores with 
radii ~ 2 — lOA and length from a few nanometers to mil- 
limeters. Biological channels Integral membrane proteins 
that form channels in biological systems are difficult to 
prepare for X-ray crystallography, accurate dimensions 
of membrane water channels are not yet available. How- 
ever, spanning cell membranes, they are ~ 50A in length. 
Based on sequential analogies with antibiotic pores and 
electron cryocrystallography |^,^ biopores such as the 
water transporting aquaporins are believed to have pore 
diameters ~ 2 — 4A. Thus, water movement inside these 
pores is statistically single-file, with "overtaking" rarely 
if at all occuring. 
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FIG. 1. Schematic of osmosis and pressure driven flow 
through membrane pores. The reservoirs (L) and (R) are 
assumed infinite. The rates ^ — {a, l3,-y, S,p, q) are condi- 
tional solvent entrance probabilities at the ends of the pore 
indicated (see text). For typical experiments, the membrane 
is impermeable to solute, AP = — P^ , (L) contains pure 
solvent, and J indicates the direction of solvent flux. 



Using a symmetric exclusion model [ p8[ , thermody- 
namics, and simple kinetic theory, we formulate an ap- 
proximate model for transport through nearly single-file 
pores which can be solved exactly. Estimates of the 
microscopic kinetic parameters used allow the theory to 
be general enough to predict trends and to explore de- 
pendences on physical parameters such as pore geome- 
try, microscopic interactions, and bulk ffuid properties. 
Our model complements, although is more flexible than 
detailed numerical simulations. Microscopic motions de- 
rived from simulations can in principle be used to evalu- 



ate the system-dependent parameters used in this work, 
provided accurate force fields are known. 

In the next Section, we review the usual linear 
phenomenological expressions used to describe osmosis 
or pressure driven flows through semipermeable mem- 
branes and identify all macroscopic variables. A one- 
dimensional symmetric exclusion model is then formu- 
lated for particle dynamics within a pore that connects 
two fluid reservoirs. The rate coefficients used are ap- 
proximated using kinetic theory near thermodynamic 
equilibrium. 

In the Results and Discussion, we explore the model 
as a function of reasonable parameters and present the 
mean flow rates in a series of plots. The macroscopic 
variables typically used are then related to the rate con- 
stants used in the microscopic theory. Nonlinearities and 
Arrenhius temperature dependences are also discussed. 
In the Conclusions, we summarize our results and intro- 
duce extensions to be studied in the following paper. The 
Appendix gives a simple result for pores with disordered 
interiors, or hopping rates. 



II. MEMBRANE AND PORE TRANSPORT 
MODELS 

A. Macroscopic Phenomenological Expressions 

Here, we briefly review the macroscopic descriptions 
which hitherto have typically been used to model osmo- 
sis. The typical starting assumption is that the solvent 
flux is a function of the solvent chemical potential differ- 
ence between two reservoirs (see Figure 1) that do not 
exchange solute. For two reservoirs with solvent at dif- 
ferent equilibrium chemical potentials. 



fc=i 



k=l 



(2.1) 



where the last equality assumes = T^, = P^, and 
Ck and dk are expansion coefficients which may depend on 
external parameters such as temperature, pore-particle 
interactions, etc. They may also depend on t he total 
solute concentration in the reservoirs. Equation (2.1) as- 
sumes linear response and expands Afj, in powers of sol- 
vent density (a virial expansion). The microscopic inter- 
actions are nearly identical for solvent in either reservoir 
(nearly all virial terms for the solvent chemical potential 
cancel) except for the infrequent encounters with solute 
in the solution (R) side. To lowest order, the (low) solute 
concentration difference gives rise to the free energy dif- 
ference which can be seen by using Xs = 1 — Xo, such 
that J can be written in terms of solute mole fraction 
difference A^s = xf ~Xs ■ ^ similar expansion in hydro- 
static pressure gives van't Hoff's result that the solvent 
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flow induced by an impermeable solute is thermodynam- 
ically balanced (J = 0) by a hydrostatic pressure differ- 
ence equal to Ap = —AH — —kBT{nlpxf' 



-ksTriTAxs, where nx 



i;p are the total par- 



ticle number densities. Small currents between (R) and 
(i) are thus assumed (fairly accurately) to be linear J « 
ip(An — Ap). Under isobaric conditions, and Xs ^ 1, 
osmotic flow is often represented by the osmotic perme- 
ability Pos (also denoted P/) where J ~ Pos{nf — n^). 

The constants of proportionality, Lp and Pos determine 
the rates of flow and have been the subject of consider- 
able attention. Measurements of membrane permeabili- 
ties and the number of pores per area, can yield single 
pore conductivities. Conversely, the number of pores in 
a cell or vesicle can also be determined if single pore 
conductances can be accurately measured or modelled. 
Fully microscopic approaches, such as molecular dynam- 
ics (MD) simulations reveal detailed microscopic infor- 
mation on how fluids particles adsorb to the pore inte- 
rior |23[2^; however, the time scales readily achievable 
by MD (^ns) cannot readily measure steady state flows 
through a long pore, or allow exploration of large param- 
eter dependent flow regimes. To achieve channel flows in 
MD simulations, an external potential (such as gravity or 
external field) is often applied to every particle [^, to 
approximate a convective velocity. However, particles do 
not experience local "pondermotive" forces (as is the case 
in osmosis or pressure driven flows) in the absence of ex- 
ternal electric or gravitational fields. In confined nearly 
ID systems, such external forces may yield qualitatively 
different behavior (such as shocks in the particle density 
1^,^) from that expected in hydraulic or osmotic flow. 

In the following, we analytically model particle dynam- 
ics through a single pore with the assumption that the 
reservoirs are in thermodynamic equilibrium, as was as- 
sumed in the expansion of J (Eqn. (2.1)) about a ther- 
modynamically defined A/iQ- We include sufhcient micro- 
scopic detail within the pore, while carefully obeying all 
thermodynamic constraints in the connected reservoirs. 



B. One-dimensional Pore Exclusion Model 

Our single-file model is similar to but more general 
than the one vacancy models of Hernandez |12| and 
Kohler [|l3|. The approach is fully microscopic, and 
does not rely on undefined hydrodynamic parameters 
(on these size scales) such as drag factors, viscosity, etc. 
Nonetheless, results from the microscopic theory can be 
associated and directly compared with those from exper- 
imental measurements. Consider the pore shown in Fig. 
which depicts a single-file or nearly single- file pore. The 
main interactions that we include are the excluded parti- 
cle (solvent) volumes within the pore; thus we adapt one- 
dimensional exclusion models |^,|l6| . The pore is divided 
into N sections labelled i, each of length £. The occupa- 
tion variable ai (t) defines whether pore section i at time 



t contains a particle {ai{t) = 1) or is empty (cri(i) = 0). 

The discrete dynamics are defined by the follow rules: 
During any infinitesimal time interval dt, randomly chose 
a particle with label < i < iV -I- 1, where i — 0, N + 1 
corresponds to choosing the (L), (R) reservoir. If i = 
(the left reservoir) is chosen during dt, and (7i{t) — 0, 
then we fill the first site (tri = 0) ^ {ai{t + dt) = 1) with 
probability adt. If i = -I- 1 (the right reservoir) and 
o'Ar(t) ~ 0, then we extract a particle from (i?) and set 
<TN{t + dt) = 1 with probability Sdt. These steps allow 
injection of the end sites of the pore from the reservoirs 
if and only if the end sites are empty at time t. If a par- 
ticle at « = 1 is chosen, we empty it into the (L) reservoir 
with probability "fdt, and if site i = 2 is empty, we move 
the particle to the right with probability pdt. Similarly, 
if a particle at site i = N is chosen, it moves into the 
(R) reservoir with probability /3c?i, and if aN-i{t) — 0, 
it moves to the left with probability qdt. If an interior 
site, 1 < i < N, is chosen, and i + l{i — 1) is empty, 
then we move the particle to the right(left) with prob- 
ability pdt (qdt) . The model assumes that p and q are 
constant along all the interior sites: Effects of site depen- 
dent p = q are derived in Appendix A. All the moves are 
allowed only if the site to be entered is empty, except for 
the reservoirs i — 0, + 1. After the a step is chosen 
and particles moved with the assigned probabilities, an- 
other particle is picked at random and the possible moves 
made. This algorithm is identical to a Monte Carlo simu- 
lation | ]30| except that spatial moves are coarse-grained to 
£ and the effective interaction energies (which determine 
the acceptance of moves) between particles are infinite 
(representing hard spheres) for a particle occupying the 
neighboring site, or zero if the site is empty. These kinetic 
steps can be summarized by the following instantaneous 
fluxes 

J,{t) = pa,{t){l - a,+i{t)) - qa,+i{t)il - a,(t)) 

l<i< N -1 



J^L)^i(t) = Mt) = a{l - ai{t)) - 7(71 (i) I = 

JN^iR){t) = JN{t) = f3<JN{t) - (5(1 - (7Ar(i)) I = N 

(2.2) 

These dynamics implicitly include effective excluded vol- 
ume interactions between solvent particles, which are 
strong in confined geometries such as ID pores. How- 
ever, we assume that solvent-solvent attractions within 
the pore are weak compared to solvent-pore particle at- 
tractions such that concerted cluster hops are rare. These 
processes can be be treated numerically |Q and is out- 
side the scope of the present study. 

To extract the steady-state cur rent s and occupations, 
we take the time average (. . .) of (2.2) to find 
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{J, (t)) = p{a, (t)) ~ q{a,+i (t)) + (q - p) {a, {t)a,+i (t)) 
(JoW)=a(l-(ai(t)))-7(ai(i)) 



o-,: = CTi - (t- Ij- = 



(2.3) 



Exact solutions to the steady state {J{t)) under cer- 
tain parameter regimes and in the thermodynamic hmit 
have been found using matrix algebra techniques |§8|,^ . 
However, in osmosis and hydraulic pressure driven flows 
across symmetric pores, where the particles do not ex- 
perience intrinsic or pondermotive forces, and have not 
developed collective stream velocities (c./. Conclusions), 
q = p. That is, a particle is as likely to hop to the 
left or right if both sites to the left and right are unoc- 
cupied. The only "force" driving osmosis is the direc- 
tionally asymmetric excluded volume along the length of 
the pore, which is ultimately determined by the entrance 
and exit rates at the boundaries i = l,N. Geometri- 
cally asymmetric pores would have p ^ q and need to 
be treated numerically. The convective stream flow limit 
where also p ^ q, is discussed in the Conclusions. Many 
molecular dynamics simulations designed to study flow 
in narrow channels impose a microscopic pondermotive 
force with p ^ q [ p5| (via a gravitational potential for 
example) in order to accelerate the flow. However, grav- 
ity and pressure forces, although macroscopically both 
derived from potentials, do not manifest themselves mi- 
croscopically in the same way. When q = p the higher or- 
der correlation {ai{t)ai+i{t)) vanishes in ( ^.3D and upon 
time averaging and invoking steady state particle con- 
servation, Ji = Ji+i = Ji_|_2 = ■ . . = Jo = Jn, where 
we have dropped the (...) notation. Since all the steady 
state currents across each interior section is identical, we 
can sum them up along the chain to obtain 



J 



P 



N - 1 



{(Ti - o-at) = Jo = Jj 



N 



(2.4) 



Equations (2.4) give three equations for the three un- 
knowns J, fJi, and (Tjv which are found to be (with 
C = {a,f3,l,S,p)) 



p{a(3 — 7(5) 



{N - l){a + -i){l3 + 6) + p{a + [3 + + 5) 

(2.5) 



and 



: : o-w = TT—c — 1 (2-6) 



a -)- 7 
which determines tr,-. 



[3 + 6 



p a + (3 

p(a(3 — 7i5) 



(a + 13) [{N - 1) (a + 7)(/3 + (5) + p(a + /? + 7 + 6)] 

{i-l){aP-l5) 

{N - \){a + -f){(3 + S) + p{a + (3 + + Sy 

(2.7) 

For an infinitely thin membrane {N — 0), the particles 
pass the mathematical surface based upon their reservoir 
kinematics alone and never interact with each other or 
the membrane interior while crossing the barrier. Over- 
taking of particles, when two adjacent, indistinguishable 
particles switch positions, does not contribute to the 
overall current. Therefore, the model above is also valid 
as long as each section i is not too wide as to contain 
more than one particle at any time. Pore diameters of 
up to 2 ~ 3 times the particle diameters will on average 
preclude multiple occupancy. 



C. Parameter Relationships 

The parameters ^ = (a, /?, 7, 5,p) used in the above ki- 
netic model are conditional probabilities and reveal the 
structure and microscopic mechanisms of osmosis and 
pressure driven flows. Here, we explicitly relate ^ to 
macroscopic thermodynamic variables. Rather than find 
precise numerical values for ^, we estimate them using 
reasonable physical models and extract their most sen- 
sitive dependences, thereby illustrating the qualitative 
features of osmosis and pressure driven flows through a 
pore. 

The main assumption we make is that the entire sys- 
tem is under local thermodynamic equilibrium (LTE): 
Locally, particle velocities obey a Maxwellian distribu- 
tion. This is justified since measurements of almost all 
microscopic systems give single pore osmotic transport 
rates of J < 10^/s. Typical pore diameters in aqueous 
solutions give mean free paths Xmfp ^ 5 A. Thus, am- 
bient thermal velocities vt — ^/ksT /m ~ 4 x lO^cm/s 
(where m is the particle mass), give mean collision times 
of Tcoii — Xmfp/vt — 0.5ps <C J~^. Therefore, particles 
within the pore equilibrate by suffering O(IO^) collisions 
while being transported from (L) to (R). This is suffi- 
cient for local thermodynamic equilibrium (LTE) to hold. 

Figures |l|(b) are schematics depicting the activa- 
tion energies determined by microscopic membrane- 
solvent interactions. Although these are predominantly 
of an enthalpic nature, determined by instantaneous in- 
termolecular potentials, a small entropic contribution 
may arise due to for example rotational averaging of non- 
spherically symmetric particles. We do not attempt to 
model in detail all the configurational dependences of 
particle trajectories, but rather assume the activation en- 
ergies to be effective quantities, averaged over say, molec- 
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ular rotations, internal vibrations, etc, which will not 
qualitatively affect the mean quantities we wish to study. 
Ep,q are internal binding energies to the pore walls which 
must be overcome when particles hop from one section 
to another and can be thought of as effective interac- 
tion energies between pore and solvent. As mentioned, 
the model presented is valid provided solvent-solvent at- 
traction within the pore ^ Ep. Ep^^ and E^^s are the 
effective energies of removing and injecting a particle. If 
we assume that there is no additional "constriction" at 
the pore ends, for a pore that repels particles relative to 
its interaction energies in the bulk reservoirs (L) and (i?) , 
Ep^^ ~ and Ea^s > (Fig ||(a)). Similarly, as shown 
in Fig ||(b), an attractive pore will have -B/3,7 > and 
Ea^s — 0. The physical origin of the energies such as 
Ea — Ery defined in Figure 2(b) are most likely from hy- 
drogen bonding interactions in aqueous systems. Trans- 
ferring a water molecule from the polar, hydrogen bond- 
ing bulk environment into the pore region would require 
breaking hydrogen bonds (~ hksT each) and forming 
van der Waals-type interactions with the pore interior. 



E(x — Efi — E^ 



L=Nl ^_ (^^- 

(b). 




-Ey -E^ 




(c). 



FIG. 2. Simplified pore-solvent interaction energy land- 
scape for (a) repulsive and (b) attractive pores. The progres- 
sion of the activation energies are indicated as pores are 
altered or conditions modified such that pore-solvent binding 
increases: (a) (b) —f (c). 

Additional activation energies at the entrance/exit 
sites would represent for example energetics of a possible 
intermediate state where hydrogen bonds are being bro- 
ken, but attractive van der Waals interactions have not 
yet formed. For simplicity of notation, we assume these 
are small and neglect them. Inclusion of these bound- 
ary activation barriers will only shift the absolute values 
of Ea, Efj, E^, Es and will not qualitatively change our 
results. 

The probabilities £, can be estimated by first passage 
time or transition state theory. The average time for 
a thermally agitated particle to climb over a barrier of 



height E^ is given by 



(2.8) 



where ^0 is a thermal attempt frequency. The transition 
probabilities per time S,dt are thus given by dt/r^. Upon 
choosing c?t -C r^. 



^~Coexp(-£;^/fcBr). 



(2.9) 



The prefactors of q = p and /3, 7, are estimated by 
POtPoiIo — vtII. Although more sophisticated transi- 
tion state theories exist, the thermal attempt frequencies 
are all on the order of Ips"^. The frequencies ao and 
(5o however, depend on the thermodynamic state of, and 
the particle numbers in the reservoirs (L) and (i?) respec- 
tively, as well as the pore entrance areas. For example, 
for effusing dense gases, ao = n^v^A^jA:, proportional 
to the number density nl^ of particles able to enter, the 
thermal velocity and the available area A^ for en- 
trance. A completely analogous expression holds for 
with (L) replaced by [R). Although the 1/4 factor is 
valid for g qualitatively similar expression holds 

for dense liquids. 

The microscopic pore area Ap can be defined by 
Eafi. First consider Ea{r) to be a general interaction 
of pushing a particle perpendicularly through a mem- 
brane. When the particle is pushed through at a posi- 
tion f away from the pore (through a lipid layer for ex- 
ample), Ea{f) ~ 00. When the trajectory runs through 
the center of an approximately circular pore entrance, 
Ea = Ea{f = 0). Therefore, we can define a pore ra- 
dius {Ap ~ Tp) by demanding activation energies such 
that trajectories outside of are energetically rare. For 
example. 



Ec,{r < rp) - E^{Q) <^lQkBT. 



(2.10) 



Condition (2.1C) simply contrains r < Tp to be within a 



region where entrance is not impeded by "wall" interac- 
tions. Thus, Ea (and Es) are the activation energies for 
a particle dragged through the approximate "centerline" 
of the pore. When more than one species occupies the 
reservoirs, but only one (the solvent) can enter the pore, 
the entrance rates will also be proportional to their mole 
fractions. For notational simplicity, we henceforth define 
a = aoXo 6xp[-i?Q/A;BT] ^"^^ ^ = ^0X0 ^^V'l-Es/ksT]. 
Thus, ao oc Ap is now the rate prefactor for pure sol- 
vent entrance (xo = 1) at a reference number density for 
a pore with entrance area Ap. Only the solvent fraction 
can enter the pore and partake in the dynamics described 
by( |2.2| ) because the solutes are too large to enter through 
Ap. 

Although q — p, a, net current will occur due to differ- 
ences in a ^ S, and/or f3 ^ j. Henceforth, for simplicity, 
we consider only structurally right cylindrical pores such 
that /3o = 7o and ao = So (or Ap = Ap and Ea = Es) 
under equilibrium (T^ = and = P^) conditions. 
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Define from the numerator of J(^, N) the dimensionless 
driving variable 



X 



a/3 — 7(5 

a/3 

l-^exp[iE^+Ep-E^-Es)/kBT] (2.11) 



where Ai; = E^-E^ + Ep-Es = E^-E^, Axs 

Xsj and Xs^^'^^' is the solute mole fractions in (i?), (i). 
The solute particles, unlike the solvents, are too large 
to enter the pore. There will be a steady state particle 
flux as long as X ^ 0, which can occur either by virtue 
of AE ^ (predominately pressure driven flow), and/or 
=xf-Xs^O (osmosis). 
Although electrostatic potential differences can con- 
tribute to AE^, they would also make q p, which we 
do not consider. Therefore, we assume that AE^ arise 
solely to a hydrostatic pressure difference in the baths 
which pushes the molecules against their repulsive inter- 
actions, and change the relative energies of activated pore 
entrance. 



III. RESULTS AND DISCUSSION 

In this section we examine the solutions of one dimen- 
sional models and discuss their physical meaning. Vari- 
ous effects of physical and chemical parameters, such as 
pore length and solute interaction dependences, are out- 
lined. 



A. Osmosis 

Consider symmetric pores (/3 = 7, E^ = Eg) connect- 
ing two reservoirs under identical hydrostatic pressures 
so that AE = 0. For simplicity, we will also assume pure 
solvent in (L) such that Axs = • Even though par- 
ticle enthalpies are identical for {L) and (i?), ^ 5q 
will drive an osmotic flow from (L) — > (i?) since in this 
limit X = (a — (S)/a ~ xf 7^ is due solely to a so- 
lute concentration difference. The steady state current 
becomes 



{N -l){a + l){a+l 



aX) + p{2a + 2-aX)' 
(3.1) 



where a = a/ (3. Since a//3 is a ratio of pore entrance to 
exit rates across the pore, it measures an effective pore- 
solvent binding constant. 

How does J(^, N) behave as pore-solvent combina- 
tions with varying a are used? When a symmetric 



pore is repelling, the energy landscape is approximately 
that shown in Fig. ||(a), where E^ ^ Es > 0, and 
Ep ~ iJ^ ~ 0. As pores with increasing solvent afflni- 
ties are used, Ea — Eg decrease while a, 5 increase, until 
the energy barrier shown in Fig. |^(b). is reached, where 
a — ao/Po. Further increasing pore-solvent attraction 
requires Ea ~ iJ^ ~ 0, and increasing Ep ~ E^. Thus, 
as pores follow the sequence (a) (b) (c), a first 
increases as exp{—Ea/kBT) (as E^ > is decreased), 
until a — ao/Po, where upon a increases according to 
eyip{Ef-j/kBT) (as i^^j > is increased). While the pore 
is repelling, and /? = /3o is constant, J(^, N) has a maxi- 
mum as a function of pore-solvent affinity at 



aiPo) 



2p/l3o + {N-1) 



-,1/2 



(7V_1)(1_X) 



(3.2) 



However, if < a*(/3o), the maximum defined by (3.2) 
is not reach by decreasing E^- Further increasing a/ f3 
requires decreasing (3 < /Sq with fixed a = ao as shown 
by the transition (b)— i>(c) in Fig. ^. The value of (3 which 
yields a maximum in J in this case is 



■f ^ - f 
" ("0) = ( ^ 



aoiN-l) 



ao{N -l){l-X)+p{2-X) 
p=0.003, X=x =0.0018 [=100mM] 



1/2 



(3.3) 
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FIG. 3. Flux as a function of effective pore-solvent binding 
a. We use representative values dt ~ lOfs, p = /?o = 0.003, 
and X = X? — 0.0018, corresponding to osmosis from pure 
solvent to a lOOmM solute solution, (a). J{^,N = 15) is 
shown for ao — 0.001(solid curve) and ao — 0.0045(dashed 
curve) which give difference maxima at q*(qo) and a*(/3o) 
respectively, (b) shows J for A'' = 5, 10, 15, 20, 25 when 
ao = 0.001. Note the discontinuities in slope at log(l/3) cor- 
responding to the point where the source of increasing affinity 
switches from the increasing of a to the decreasing of /3. (Fig- 
ure 1(b).) 

Figure |^ shows the current as the solvent-pore affinity 
is increased according to a, as a = ao//? ^ ao//3o(= 
ao) ~^ ci/Po- The solute concentration is represented 
in molar units [cs] = xf ^ 55.56, such that a lOOmM 
aqueous solute solution corresponds toX — xf — 0.0018. 
Although J and ^ can all be rescaled by the constant p, 
for concreteness and comparison with experiments, we 
consider a time step dt ~ lOfs <JC Tp « 0.3ps, which 
sets /3o ~ p — 0.003 and yields flows in the range of 
10* X J (particles per /is) in good agreement with Pos ~ 
10~^'*cm^/s measured across aquaporin water channels 
with physiological lOOmM) osmolyte solutions [HJ]. 

Figure ^(a) plots the flux (/xs~^) as a function of 
log a for two different values, ao/Po = 0.001/0.003 < 
(a//3o)* = ^8/7(1 -X) and ao//3o = 0.0045/0.003 > 
(a//?o)*. At log(l/3), the curve bifurcates into the two 
solutions corresponding to either de crea sing a or increas- 
ing /3, with maxima determined by (3^) and ( |3.3[ ) respec- 
tively. The discontinuity in slope at log(l/3) (dashed 
curve) is apparent, although the one at log(3/2) near 
the maximum is not. Figure |^(b) compares J(^, N) for 
vario us pore lengths L — N£ at ao — 0.001 (such that 
( |3.3D pertains). Whether ao//3o is greater or less than 
a*(/3o) depends on the reservoir number density and Ap; 
for gase s, ao <^ /3o a-nd a curve with maximum defined 
by ( |3.3| ) obtains. However, dense liquids entering large 
pores may have relatively large entrance rates such that 
(^o/f^Q ^ a*(/ 3o) and J (£, N ) may have a maximum de- 
fined by (|3.2|) . From ( |3.2|) , this latter condition is less 
likely for larger driving X, and thus a*(/3o). 



2.0 4.0 6.0 8.0 10.0 




(a). 



o/X=0.0018,0.18), a„('X=0.0018) 




0.00 

2.0 4.0 6.0 8.0 10.0 
a/p [=affinity] 

FIG. 4. Occupations a for p = /3o = Qo = 0.003. (a) 
Pore occupations for X — 0.0018 are all indistinguishable on 
this scale (solid curve). Only at X = 0.18(10M solution) 
is (jjv=3o (dashed curve) sufficiently different from ai. (b) 
100x(CTi - o-jv) as a function of a/ (3 for p = 0.003,0.03,0.3. 
Note the different behavior when a = ao = 0.003 at 
a/P — 1.0. The dashed fragments represent the occupation 
differences if ao > 0.003. 

Figure ^(a) shows the average occupation ai{a) for 
X = Xs — 0.0018. The dependence on spatial position i 
is weak and is not apparent on the scale of Fig. |4|(a), ex- 
cept if X = 0.18 (corresponding to lOM solutions) when 
(Ttv can be distinguished from cti. Therefore, from Eqn. 
( p.6D a/{a + f3) is an approximate value of averaged oc- 
cupations at small X. The small difference in occupa- 
tion along the chain is consistent with the assumption 
of local thermodynamic equilibrium. Figure ^ shows the 
difference 10'^ x (cri(a) — cr7v(a)) with ao/(3o = 1.0 (solid 
curves). The dashed curves represent the occupation dif- 
ference 10^ X (di — ctat) if ao/Po 3> 1. 

The behavior of J(^, N) as a function of affinity a can 
be further understood in terms of the associated pore 
occupations. At low a, the pore is repelling and con- 
tains few if any particles that can be transported. At 
high affinities, and occupations, most internal hopping 
steps are not fulfilled due to occupied neighboring sites 
choking off t he flow. Only at intermediate a and near 
( |3.2[ ) or ( |3.3[ ), and intermediate occupations will maxi- 
mum flow occur. At larger X, decreases (see Fig. 
^(a)), decreasing overall pore occupation and increasing 
the critical afhnities a*(/3o) and a*(ao) where the max- 
imum occurs. Therefore, the condition represented by 
.3[) becomes more relevant for larger osmotic pressures. 
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dt=10fs, p=^g=0.003, a^=0.002 






0.0 0.2 0.4 0.6 0.8 1.0 5 10 

X=x, (Molar) 

FIG. 5. The driving force dependence of J(^, 20) for osmo- 
sis from pure solvent to solution with p = /3o = 0.003 and 
Qo = 0.002. Various values of a = a//3 are shown, (a) For 
solute concentrations < 5M, the current is linear in X = x.?- 
(b) At higher solute concentrations, nonlinearities in J(^, 20) 
arise, particularly for high values of a an pore occupation. 

Figures || show the particle flux per microsecond, 
10^ X J{^,N = 20), as a function of X ~ xf at 
Po = P = 0.003, ao = 0.002 for various pore-solvent 
affinities a. For the low values a = 0.01(dashed curve) 
and a — 0.1 (dotted curve) the pore is repelling, while for 
a = 1.0(solid curve) and a = 100(thick curve) a — aQ 
while P decreases from (3q. The solid curve with a = 1.0 
is near the maximum defined by (3.3) while a = 100 is 



within the choked flow limit represented by the decreas- 
ing branch on the high a side of Figures |(a,b). The flow 
is linear for low X as shown in ||(a) . Figure g(b) clearly 
shows that nonlinearities are important for large a and 
that the flow for a — 1.0(solid curve) is lower that that 
for a = 0.01(dashed curve) only for small driving X. As 
the nonlinearity sets in for a = 1.0, the flux increases 
faster than that for a = 0.01 and surpasses it near X^^ 3. 
The flow nonlinearities shown in Fig |^(b) arise at large 
solute concentration differences and are only those result- 
i ng f rom the nonlinear dynamics determined by equations 
(2.2) and do not include solute nonidealities, unstirred 
layers, etc. That is, the coefficients 



a''p[{a + l)iN -l)+p] 



k-l 



(3.4) 



{a + 1)'' [{a + 1){N - 1) + 2pf 

contain only the kinetic parameters describing the 
solvent-pore interactions. Furthermore, our treatment 
for a and /3 assume thermodynamic equilibrium, or no 
unstirred/polarization layers. The modifications neces- 
sary to include unstirred layers are treated in the follow- 
ing paper. 



From the denominator of (3.1), J is linear in X for 
aX/{a -I- 1) < 1. From (|2!|), we see that a/{a + 1) is 
the average occupation for small J. Therefore, nonlin- 
earities are important only when occupation at is high, 
when exclusion interactions, and thus nonidealities, are 
most pronounced (the a — 1 .0, 100 curves in Fig. ||(b)). 
Although Eqn. (pTl|), or (3.4) fiflly describes particle ex- 
clusion nonidealities within the pore, for the small Axs 
usually encountered, a linear relationship. 



P^Xs 



d + l{N - l){a + l) + 2p 



(3.5) 



is accurate as demonstrated by Fig ||(a). In the limit 
p ^ {a + 1)N, Pos approaches 



Pn. 



ap 



[a + l)^nT' 



(3.6) 



In this limit ( |3.q ), the "rate limiting steps" are the 
internal steps (p) of pore particles, and the overall 
flux scales as In the limit represented by ([3.(^), 

an Arrhenius temperature dependence of {Ea — Ep — 
Eq)/kBT[{Efi - Ea - Eq)/kBT] is expected for a > 
l[a<C 1]. 

" (X ap (X 



According to the form of ao i Pos cc ap (x for a ^ 1 , 
r^'^ for a ^ I. Note the curious 



and Pos (X p/a oc 
Pos fx dependence indicating the a > a* regime 
shown by the large logj^o ^ regions of Fig. ^ where de- 
creasing Tp (and a) increases J by decreasing the high 
pore occupation. 

Now consider the limit p 3> (Q;-|-l)iV. The rate limiting 
steps now involve pore entrance or exit and 



Pr. 



a 



2{a + l)nT 



(3.7) 



independent of L. This limit is equivalent to the single 
site {N — 1) model since only the boundary entrance and 
exit rates are relevant. The radius and temperature de- 
pendences in this case are Pos oc a oc exp(— E'^/fesT) 
and Pos (X f3 (X exp{—E0/kBT) for a <^C 1 and a ^ 1 
respectively. Thus, the pore length L ~ N£ can deter- 
mine the temperature dependence of J since it deter- 
mines which of the two limits (3.6) or (3.7) is relevant. 



Except for experiments where lipids undergo phase tran- 
sitions which affect the pore structure |Q , a single expo- 
nential temperature dependence is almost always experi- 
mentally observed in biological systems involving osmosis 
p2t , indicating one of the above limits for Pos holds in 
the regimes explored. One interpretation is that a num- 
ber of hydrogen bonds (~ bhsT each) must be broken 
before water can enter the pore. This is consistent with 
the identiflcation of Ea as the ener gy re quired to break 
water H-bonds in the a <C 1 limit of (3.6). Fluxes of non- 



polar solvents through artiflcial pores may be expected 
to display richer temperature dependences since a wider 
variety of temperatures and molecular interactions can 
be accessed and experimentally probed. 
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B. Pressure Driven Flow 



where in this case, 



For incompressible reservoirs without solutes, Axs = X — 1 cxp [{Ea{P^) + Ep E~^{P^) Es)/kBT] 



Xs — xf — 0; however, flow can be driven by differences 
in hydrostatic pressure. Hydrostatic pressure variations, 
rather than controlling the rates via the number density 
in the pre-exponential factors ao, (Sqj modify the relative 
activation energies For example, consider increas- 
ing keeping P^ constant. For nearly incompressible 
liquids, the particle density changes negligibly, but their 
particle enthalpies increase as they are pushed against 
their interaction potentials. Changing the relative reser- 
voir and pore enthalpies changes the entrance and exit 
kinetic coefficients according to Fig. ^(a), (b) for initially 
repelling or attractive pores. 



,L E„=0 



Ey(P^)>0 



E.=0 



(a). 



pL E„=0, (P^- )>0 



E.=0,E,>0 



Ey=0, E^ {pL )>0 



(b). 



FIG. 6. Hydrostatic pressure increases in left reservoir, (a) 
In attractive pore, a — 5, Ejs is constant, and i?^ = Erj[P^). 
(b) For repelling pore, and small pressure differences (lower 
solid and dotted levels) 'y = f3, Es is fixed, and Ea = Ea{P^). 
Pressure differences larger than Eg result in Ea « and 
E-, = E^{P^) > 0. 

Since we assume only P^ is varied, Es, Ep are kept con- 
stant, while Ea = Ea{P^) and E^ = E^{P^)- First 
consider an attractive pore where only E^{P'") and 7 
changes as P^ is increased, as schematically depicted in 
Fig. p|(a). The flux is 



Jh{P^),N) = 



{N - l){a + 1 - X)(a -h l)+p(2a 2 - X) 



where 



X(F^) = 1 - exp [AE/kBT] 



(3.8) 



(3.9) 



and AE = Ef3- E^{P^). 

When hydraulic pressure pushes particles through a 
repelling pore, the current can be written in the form 



j(a(P^),7(P^),A^) = 



p(5X 



{N - \){5 + \){5 + 1 - X) +p (7 + ^ + + 1)(1 - X)/7) 

(3.10) 



(3.11) 

and d measures the solvent-pore affinity since a is varying 
with pressure. For low pressure differences such that 



pL _pR ^ pj^ 



fdAE\ 



(3.12) 



flow across repelling pores can be further simplified since 
7 = /? as evident in the lower three pressures P^ in Fig. 
^(b). When the pressure pushes the enthalpy of the left 
bath above that of the pore, the full pressure dependence 
( |3.10D is required. Equations ( ^ ) and (3.10), along 



with the equation of state of the solvent in the reser- 
voirs completely determine the nonlinear flux-pressure 
relationship. However, for low pressures, a Taylor ex- 
pansion about — P^ yields 



X « - 



AP l'dAE\ 
\dAp) 



T,AP=0 



. AP 



(3.13) 



where the last equality is valid for incompressible liquids 
and V is the molecular volume of a solvent particle in the 



reser voirs (for H2O, 5 ~ 3 x 10 ^cm'^). Linearizing (|3.8|) 



and (3.10) about small pressure differences, and recalling 
that a ~ i5, we find 



Lr. 



p 

^ OS 



(3.14) 



as required by linear noncquilibrium thermodynamics. 
However, differences between hydraulically driven and 
osmotically driven transpo rt a ris e a t highe r ord er in X, 
as expficitly calculated in (^.5]), (^), and ( 3.1C| ). 

Continuum Poiseuille flow (zero Reynolds number so- 
lution of the Navier-Stokes equation for flow through an 
infinite right circular cylinder |]l9| ) expressions (where 
J oc Tp) have been proposed to describe transport 
through microscopic channels. Although these contin- 
uum relations have no validity in microscopic settings, 
they are nevertheless often used to obtain estimates of 
permeabilities, pore radii pO[ |, and membrane pore den- 
sities. These estimates typically do not agree well with 
independent measurements of pore radii and pore perme- 
abilities . Although the Tp dependences in our model 
are valid only for a limited range of pore radii, they are 
consistent with the inapplicability, at microscopic scales, 
of the continuum expression Pos oc r'^/{rjL), where 77 is 
the bulk solvent viscosity. 

Fluid flow through small orifices has been described 
as being both "viscous" and "diffusive." Even though 
the limit (3.6) appears to be "viscous" due to the L^^ 
dependence, and and (^) "diffusive," we see that both 
pressure and osmosis driven fiows in the present model 
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is diffusive in nature. Consider the time averaged rate 
equation for occupation at site i, 



= -p(Ti{l - CTi+i) 
+P{1- 



(3.15) 



which on scales » ^ is a diffusion equation with pf^ = 
a cooperative diffusion coefficient. When particles in- 
teract along the chain the overall flux is determined by 
the cooperative diffusion of particles. In the simple one- 
dimensional model presented, cooperative diffusion is re- 
lated directly to the precisely defined dynamical rules 
( ^.2D , but in general is defined by linear density-density 
correlation functions, or Green-Kubo relations. 

The neglect of true convective, or momentum transfer 
effects can be neglected in the systems we consider. The 
signature of convective flow, such as Poiseuille flow, is a 
particle velocity distribution function with is Maxwellian 
centered about a small finite stream velocity. However, 
in ID, or nearly ID pores, the particles thermally equi- 
librate with the stationary pore walls very quickly and 
are unable to collectively transfer a stream momentum. 
This is particularly valid for repelling pores where en- 
tropic pressure results in low pore occupations cr^ ^ 1. 
In this large Knudsen number regime H, the particles 
will collide with walls more often that with neighboring 
particles. 



IV. EXTENSIONS AND CONCLUSIONS 



An extension of the ID chain represented by (2.2) to 
include a distribution of site-dependent internal hopping 
rates pi — qi can be straightforwardly calculated, as out- 
lined in the Appendix. Biological realizations of inter- 
nal pore defects with site dependent hopping rates are 
gramicidin A channels comprising of two barrel struc- 
tures joined at the defect (membrane bilayer niidplane) 
101 . Similarly, water permeation through a bilayer mem- 
brane can be viewed as pores along the hydrophobic 
lipid tails with the bilayer midplane and possible un- 
saturated bonds acting like defects. Channel proteins 
with widely varying internal interactions would also be 
modified according to the results found in the Appendix. 
Strong particle-pore interactions which are incommensu- 
rate with i may also lead to nonlinear transport. Here, 
possible cluster diffusion |^l| and locked-sliding phases 
such as that found in Frenkcl-Kontorowa |Q and fric- 
tion models |35j can also affect transport. These effects 
offer the possibility of interesting nonlinearities amenable 
to numerical study. Solvents that have strong mutual at- 
tractive interactions can also be considered. For example, 
a particle may hop with different rates into an empty site 
at i -|- 1 depending on whether or not there is an attrac- 
tive particle behind it, at site i — 1. The probability of 
such a hop to the right from i ^ i + 1 can be represented 



by ptTi(l - (7i+i)(l - (Ti_i) +p'ai{l - a^+i)ai-i where p 
represents the rate of an isolated particle, while p' p 
represents the rate in the presence of another particle at 
i — 1 binding to and preventing the hop of the particle 
at i to i + 1. The resulting particle density profile along 
the chain will no longer be linear. These effects offer 
the possibility of interesting nonlinearities amenable to 
numerical study. 

For pores that are slightly larger, with radii larger 
than a few times the particle diameters, similar stochas- 
tic treatments can be applied provided particle momenta 
relax much faster than particle positional degrees of free- 
dom. Only then is the collective stream velocity zero, and 
p — q can be invoked. Near the single file limit, the pore 
walls are in constant interaction with the transported 
molecules and randomizing their momentum. Neverthe- 
less, for larger pores, we can estimate the relative rates 
of positional and momentum relaxation. Pressure (ei- 
ther hydrostatic or osmotic) driven fiow through a macro- 
scopic circular pipe is described by Poiseuille's law p9[ 



(4.1) 



where v{r) is the fluid element velocity along the pipe, 
and T] is the bulk fluid viscosity. Equation has of- 
ten been extended to describe flow through microscopic 
pipes, especially biological membrane pores, having radii 
Tp of Angstroms. This is the origin of the notion that 
osmosis through microscopic pores is "convective." How- 
ever, for the Poiseuille stream convection to be compara- 
ble to "driven diffusive transport," the ratio 



Pe^^^- 



p 



1, 



(4.2) 



where Dc is the effective collective diffusion coefficient 
which describes dissipative particle density relaxation, 
and is also approximately the momentum diffusion, or 
sound attenuation coefficient. For a hydrostatic pressure 
difference of AP — 2.5 atm (corresponding to a lOOmM 
solute concentration difference), and bulk viscosity and 
attenuation coefficient 77 — O.Olg cm/s, and Dc — r]/p, 
convective, or "viscous" flow will be irrelevant for 
Tp ^ 1000 A. For typical zeolites and biological pores 
under the above conditions, r„ ~ 3 — 8 A, Pe ~ 10~^ <C 1. 



Thus, "convective" flow described by (4.1) is negligible 



compared to diffusive transport in systems of present con- 
cern. Diffusive transport in this context should not be 
confused with tracer diffusion or collective diffusion in 
binary mixtures, rather, it is the consequence of micro- 
scopic momentum relaxation due to frequent collisions 
with the pore interior. These collisions increase relative 
to particle-particle interactions for smaller and smaller 
pores, where the inner surface to pore volume ratio in- 
creases. 

Additional support for the basic assumption of diffu- 
sive particle motions is found from nonequilibrium molec- 
ular dynamics studies (NEMD) of flow through a slit |lj] 
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where an external driving force q ^ p was not used explic- 
itly. However, the two reservoirs were held at different 
chemical potentials, or pressures. The density profile for 
Lennard- Jones particles with attraction energies close to 
/cfiT, driven across SOAlong slits of widths 2.5 times the 
particle diameters with a 35 atm hydrostatic pressure 
diffe rence was found to be linear, suggesting validit y o f 
( ^.6| ) and Pick's law even under extreme conditions [Q . 
The validity of (4.2) and diffusive {p = q) flow suggest 
that pores with radii '--^ 10 — 50A a few times molecu- 
lar diameters can also be treated by multidimensional 
exclusion models consisting of a series of ID chains with 
particle interchange among them. Solvent-solvent attrac- 
tive interactions and diffusion of small clusters can also 
be straightforwardly implemented. 

Although we have treated only a single pore, actual 
experiments measure flow over an area containing an en- 
semble of pores. The averaged steady state flux neglect- 
ing unstirred layers would then be 



(J)= /(7V,{p.})J(iV,fe}) 
Jv,{p.} 



(4.3) 



where f{N, {pi\) is the number distribution for pores of 
length iV, with the set {pi} of internal defects. The 
pores in ( [1.3| ) need not be straight, and the sum over 
N ~ d/ ^ oo represents a distribution of pores with 
various total arc-lengths traversing a membrane of thick- 
ness d. We have not considered the nonlinear effects of 
solute-pore interactions, which we treat in detail in the 
companion paper. 

The main results of this study is a thermodynami- 
cally consistent model of osmosis and pressure driven 
flows through nanoscopic-sized pores. All parameters in 
the theory are microscopically motivated yet are general 
enough to allow predictions in trends and identify cor- 
relations between transport rates and macroscopic ther- 
modynamic properties. Predictions from this transport 
model are the site occupancies (2.6), the dependences of 
X on Axs and AP = — P^, the Arrhenius temper- 
ature dependences, the por e radius rp dependences, and 



the flow maximum (3.2-3.3) as a function of pore-solvent 
binding constant d. Also, microscopic expressions are 
derived for phenomenologically measured permeabilities 
Pos and Lp. Equation (4.S), and the microscopically 
estimated values for the parameters involved offer an 
reasonable, semi-quantitative description of osmosis and 
pressure-driven flow through microscopic channels. The 
analytic model presented provides a simple, physically 
motivated analysis that complements more complicated 
numerical simulations by predicting rich parameter de- 
pendences. 
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APPENDIX A: FLOW THROUGH DISORDERED 
PORES 

The exact steady state flow through a pore with ran- 
dom internal hopping rates Qi = pi can be found by 
extending the uniform pore solution (3.1). When these 
random hopping rates are time independent, or fluctuate 
independently of the particle occupations a, the steady 
state between section i and i -I- 1 is 



(Al) 



and the mean-field solution becomes exact. Summing 
( |Al| ) across the channel yields 



N-l 



which upon solving with (2.6) for J, 
a/3 — 7(5 



J 



{a + 7)(/3 + 6) EtiPT' + {a + f3 + -/ + S) 



(A2) 



(A3) 



The above result also constitutes a steady state solu- 
tion to the XXZ spin chain Hamiltonian with site- 
dependent energies. 
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